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1. Introduction 

Rational surfaces play an important role in the frame of practical applications, 
especially in Computer Aided Geometric Design (see ^1 El and the references 
therein). These surfaces can be parameterized, i.e. can be seen as the image of a 
generically finite rational map 

(j): K2 K3 

where K denotes the ambient field, which we assume to be of characteristic zero. 
In the sequel we will address the following questions: 

• Properness problem: decide if (f>, co-restricted to its image, is invertible. 

• Inversion problem: If (f>, co-restricted to its image, is invertible, then com- 
pute its inverse. 

Both questions have already been solved theoretically and algorithmically. For 
plane curves, the situation is very well-known: one can relate the properness and 
inversion problems to Liiroth's theorem, and there are different algorithmic pro- 
cedures to solve them (see [131 [17. .21. .22) '). For space surfaces, there exist some 
algorithmic approaches based on u-resultants [7j and on Grobner bases in In 
|18| a complete algorithm is given to solve both problems by means of univariate 
resultants and GCD computations. Our starting point is the resultant matrix-based 
method presented in [131 Chapter 15] for inverting a parametrized algebraic sur- 
face, and in |31 §5] where it is used for computing the inverse image of a point of a 
parameterization. 

In this paper, we introduce a general matrix-based approach for dealing with 
both problems. We will begin by reviewing the plane curves case; this will clarify 
our approach and help the reader to understand it. In section |31 we will introduce 
inversion matrices associated to parameterizations. We will show that their exis- 
tence implies the properness of a given parameterization and also that they can be 
used in order to produce an inverse map in terms of quotient of determinants of 
some of its sub-matrices. 
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The problem of deciding whether a given parameterization is proper or not is 
not an easy task, hence in general the construction of such matrices is a non trivial 
problem. In ^ some matrices called Jacobian matrices were introduced in order to 
deal with the properness problem. This is a first example of a family of inversion 
matrices; we will briefly recall this construction at the end of section O 

In order to obtain other groups of inversion matrices, wc introduce implicitization 
matrices in sectional They are essentially square matrices whose determinant is 
a nontrivial multiple of the implicit equation of the given surface (although we do 
not need to compute the implicit equation in the process). We will prove that 
they characterize the properness of a given parameterization and moreover produce 
inversion matrices if the parameterization is proper. Then, we give a list of known 
algorithms to construct implicitization matrices. We end the paper with some 
illustrative examples. 

We would like to emphasize that the scope of this paper is not to provide a 
complete algorithm to solve both properness and inversion problems, but to give 
general theorems showing that it is possible to test properness and to get inverse 
maps of a parameterization by using tools from elimination theory which have 
been developed in principle for other purposes: as soon as one has what we call 
an inversion matrix (or an implicitization matrix) it is possible to read almost 
immediately the inverse map from it. 

A complete algorithm for computing the inverse map from the parameterization 
by using our approach requires only the construction of an implicitization matrix. 
So, the description of such an algorithm is only the description of the construction 
of an implicitization matrix. This is why we do not deal with algorithms in this 
paper, they can be found in the literature of matrices from elimination theory. 

In section 0] we give a list of some general constructions yielding implicitization 
matrices and it turns out that this list already covers a lot of cases. For the general 
situation, no process is yet known (as far as we know), but one can also try to 
produce it case by case by computing particular syzygies. 

It should be pointed out that our approach is very different from what has 
been done in 13 1181 since we provide inversion formulas in terms of quotient 
of determinants of sub-matrices of a given matrix instead of producing expanded 
rational symbolic expressions. For instance, as a consequence of our results, we 
will show that, in the case where resultant matrices give implicitization matrices 
(this is always the case for plane curves for instance), the inversion map can be 
represented by sub-matrices of a resultant matrix which is built from the coefficients 
of the given parameterization; in other words, no symbolic computations are needed 
if wc stay at this level of representation (note that this kind of representation in 
Computer Aided Geometric Design as already been widely used, for instance to deal 
with the surface/curve or surface/surface intersection problems). Moreover, it is 
interesting to know that certain big polynomials can be represented as determinants 
of certain matrices. In this paper we prove that inversion maps are given in terms 
of determinants (as soon as one has an implicitization map). 

Finally, it should be mentioned that the two main theorems 13 . 21 and 14 . 21 proven in 
this paper are also valid for rational parameterizations of hypcrsurfaces (the proofs 
work verbatim), that is to say for rational parameterizations from K"~^ to K", 
with n > 2, whose closed image is a hypersurface. We chose to stay in the context 
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of surfaces because of their important applications in Computer Aided Geometric 
Design. 

2. Preliminary: inversion of parameterized plane algebraic curves 

Before dealing with the case ^ of surfaces wc first quickly describe our approach 
to properness and inversion problems in the case of plane algebraic curves for the 
sake of clarity and completeness. 

Suppose given a rational plane algebraic curve C parameterized by 

where we can assume, without loss of generality, that gcd(pi, qi) = gcd(p2, 92) = 1- 
We moreover assume that m := max(deg(pi), deg(qi)) > 1 and similarly that n := 
max(deg(p2)j deg(g2)) > 1; if it is not the case, then C is a line and the properness 
and inversion problems are easy. 
It is well-known that 

Res,nAPi{t) - ^qi{t),P2{t) - yq2{t)) = C(x,y)'^°s(^) 

where C{x,y) denotes an implicit equation of the curve C and ReSm,n(^,^) the 
classical Sylvester resultant. This resultant can be computed as the determinant 
of the so-called Sylvester matrix Sm.n{pi{t) — xqi{t),p2{t) — yq2{t)) which satisfies 
the equality 



*S'm,n(pi(t) - Xqi{t),p2{t) - 2/(72(0) 



/ t^-\pi{t) - xqi{t)) \ 

t{piit) - xqi{t)) 
pi(t) - xqi{t) 
t"'-\p2{t) - yq2{t)) 

t{p2{t) - yq2{t)) 
V P2{t)~yq2{t) ) 



We denote by M the sub-matrix of the above Sylvester matrix obtained by erasing 
its last column. For alH = 1, . . . , m+n, we also denote by the signed determinant 
of M obtained by erasing the i^^ row. In this way, we have 



(2) 



Resm,n{pi{t) ~ xqi{t),p2{t) - yq2it)) = ^ qA^, 



where the c^'s are the entries of the last column of the Sylvester matrix (and are 
hence polynomials in IK[y]), i.e. P2{t) — yq2{t) = X^S"" Cit"^^"^^^\ 

Proposition 2.1 (^ section 2]). With the above notation we have 

deg{(j)) = 1 ^ gcd(Ai, . . . , A,„+„) e K \ {0}. 

Moreover, if dcg(t/)) ~ 1 then for all i ~ I, . . . ,m + n — 1 the rational map 



A,_^i 



is an inversion of (p. 
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It is important to notice that this proposition gives closed universal inversion 
formulas, that is to say that the inversion maps can be pre-computed in terms of 
the coefBcicnt of the polynomials pi, qi,p2 and (72- Moreover, if we stay at the level 
of matrices, no symbolic computations are needed: the inversion formula is just a 
quotient of two determinants of sub-matrices of the Sylvester matrix which is itself 
directly built from the coefficients of the polynomials pi, qi,p2 and 92- In this way, 
we can represent an inversion map by two sub-matrices of a Sylvester matrix; when 
an inverse image is required we just have to instantiate the corresponding point in 
these sub-matrices and then compute the quotient of two determinants of numeric 
matrices. This is very similar to the fact that the Sylvester matrix represents the 
implicit equation of C. 

Finally, note that we could use in the above process the Bezout matrix instead of 
the Sylvester matrix that we chosen for simplicity. For a definition and properties 
of the Bezout matrix, see 

Example 2.2. Consider the following easy example of the unitary circle parame- 
terized by 

2t 1 - f2 ^ 



1 1 



1 + t2 ' 1 + t2 



The associated Sylvester matrix is 



from we extract the matrix 







-X 





2/(1 + 




2 - 

— X 


-X 

2 






{ - 


~x 




- 


-i-y\ 




2 


— X 







— X 


2 






\ 


— X 


) 





The 3x3 minors of the matrix M are then 

Ai 2x{y ~ 1), A2 = -2x2, ^ _2a;(?/ + 1), A4 = 2x^ - 4{y + 1). 

Their gcd is a constant so we deduce that 4> is proper and we can check that the all 
the inversion formulas given in proposition are equal, that is to say that 

l\2 ^3 ^4 

where Iq = {x"^ + y^ — 1) is the defining ideal of C; for instance: 



Ai 
A2 



A2 
A3 



2x{y-l) 



-2x^ 



1 



= e Yva.c{K[x, y]/Ic). 



-2x^ -2x{y+l) x{y + l^ 

We can also check we got inversion formulas, that is to say that 

VKAi) _ VKA2) _ VKA3) _ , ^ ^f.. 

^(A2) ~ ^(Aa) - ^(A4) ~ ^^'> 
where ip is the map K[x,?/] K(<) : P[x,y) >-> P(^^-||y, ^||jy)- For instance: 
■(/-(Ai) fl + t^\ ( , \-t^ 



V'(A2) 



2t 
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In the following, we extend this approach to the more intricate case of algebraic 
space surfaces. To this task, we will introduce two distinct notions: 
Inversion matrices: these matrices will play the role of the matrix M above. As soon 
as such a matrix exists we will prove that the given parameterization is birational 
and we will deduce inversion maps similarly to proposition |5J 

Implicitization matrices: these matrices will play the role of the Sylvester matrix 
above. From such a matrix we can characterize the properness of the given pa- 
rameterization (j) and then produce an inversion matrix which is then used to get 
inversion maps. 

3. Inversion by means of matrices 

Suppose given a rational surface S parameterized by the map (0). In this sec- 
tion we develop a matrix-based approach to the inversion problem. To do this we 
introduce a certain class of matrices associated to parameterizations that we will 
call inversion matrices. We prove that if such a matrix exists for a given parame- 
terization, then this parameterization is birational and we can derive an inversion 
map in terms of sub-matrices of this inversion matrix. 

From now on, we will turn to projective geometry, so we will assume that our 
given rational surface S is embedded in P'^ and parameterized by 

0: P^ 
^' t := (ti : t2 : is) ^ (pi(t) : p2(t) : PsW : P4(t)) 

where the pi(t)'s are homogeneous polynomials of the same degree. Also, we 
will hereafter denote by K{S) (resp. ]K(P^)) the field of rational functions over 
S (resp. over P^) (see e.g. [HI §1-3.]). 

We recall that a moving surface of bi-degree (to; n) is a bi-homogeneous polyno- 
mial in the sets of homogeneous variables t := (ti,i2,i3) and X := (Xi, X2, X^, X4) 
of degree to > 1 and n > 1 respectively. Thus a moving surface M{t;2£.) can be 
written as: 

M(t;X) = Yl Mt)X", 

\a\—n 

where a G N"* is a multi- index and Aa{t) is a homogeneous polynomial of degree 
TO in t. We will say that this moving surface follows the parameterization if M 
is identically zero after substituting Xi by Pi{t) for i ~ 1,2, 3, 4. 

Definition 3.1. A matrix M is called an inversion matrix of if it satisfies the 
two following conditions : 

(i) There exists an integer to > 1 and a subset V = {t"i, . . . ,t"''} of the 



Ah{t;X) 



Md-i{t;X) 

where the polynomials A/,;(t; X) are moving surfaces following (f> of bi-degree 
{m;di). Moreover there exists (t^^ , t'^^ t'^-^ ) € satisfying tit* = tsf^i 
and tsf^-^ = i^^l}2 f^g^nd hence fit* = fa^^^ ). 
(ii) The rank o/M over the field K(S') is exactly d ~ 1. 



monomial basis {t", 




\a\ = 




/ t"i \ 






t"^ 




*M 








V ) 
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An inversion matrix M is thus a non-square matrix of size d x (d — 1), where 
d > 3. Note that the condition (ii) they have to satisfy can actuahy be checked 
exphcitly; one way to do this is described in ^ proposition 3.5 and remark 3.6]. 

Before stating the main result of this section we need to recall carefully what 
is the degree of the map 0, map that we now assume to be co-restricted to S. 
Roughly speaking the degree of is the finite number of preimages by </> of a 
sufhciently generic point on S. More precisely, the map induces an injective 
morphism 0" : K{S) -> K(p2) : / -> / o 0. Thus IK(p2) is a finite extension field 
of ]K(iS') and its degree is, by definition, the degree of (/>; usually this is summarized 
by the formula deg((/)) := []K(P^) : IK(S')]. It follows that is birational if and only 
if (f>^ is an isomorphism, i.e. has degree 1 (see e.g. j2I §1-4.]). 

Theoreom 3.2. With the above notation, if there exists an inversion matrix M of 
(j) then (j) is birational. Moreover, denoting by Aq, the signed minor of M obtained 
by erasing the line indexed by the monomial t", an inversion of (j) is given by the 
rational map 

5 ^ p3 ^ p2 . ^ ^ (A,3, {X) : A^, [X) : Ap^{X)), 
where /3i,/?2 o,nd (5^ are such that tit^^ = t^t^^ and t2tP'^ = t^t^'^ . 

over IK(S') is d — 1 wc deduce that the 



Proof. Since we assumed that the rank of 
kernel of *M is generated by the vector 



(4) 



An 



An 



An 



One has, for alH = 1, . . . , d, 
signed {d—1) x (d — 1) minor of 
But by definition we have, in K[t] 

(5) V"(M) : 

Consequently both vectors 

( t"i •■• t"'' ) and ( 0«(A„,(X)) 



Aq,;(X)) ~ Aq. ((/)(t)) and hence it equals the 
obtained by erasing the line indexed by t"'. 



/ A/i(t;0(t)) 
V Mrf_i(t;0(t)) 



/ 



V 



(Ao,(X)) 



generate the kernel of * 
by an invertible element of 



over ]K(P^), and thus equal up to the multiplication 
^). In particular we have 



t_i 

t2'' 



which imply that 0" : K{S) 

Now let the rational map 
an inverse of (p and -0" : K(P^) - 
of M we have in K[S'], and hence in K(5), 

/ 



and 



t_2 

tz' 



Aft®; h 

^) is an isomorphism. 
5- C P3 ^ P2 : X [ipiiX) : '4>2{X) : i'z^X)) be 
K(S') its associated field embedding. By definition 



(6) 




MM{X)-X) 



But for all i we have 

0«(M,(0(X);X)) = 0«(Af,(^«(t);^))) = Af,(0« o^«(t),0»(X)) = A/.(t,0(t)) = 
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in K(P^) and hence Mi{'ip{X_); 20 = in K(S'). Therefore we deduce that the vector 

(7) ( ij{x)"^ V-®"" • • • ^l^ix)°"' ) 

generates the kernel of *M* over K{S) . On the other hand we know that the vector 
Q is also a generator of the kernel of *M* over K(S'). It follows that vectors Q 
and Q must equal up to the multiplication by an invertible element of K(S'), and 
the claimed result is proved. □ 

In general, the computation of an inversion matrix is not obvious. We will mainly 
obtain them from matrices coming from elimination theory that we will describe 
in section 2J There we will sec that in most of the cases there exist algorithms 
for constructing inversion matrices. Moreover, we will also see that we can take 
advantage of the matrix formulation we have for producing closed and universal 
inversion formulas for some classes of surfaces by using resultant-based matrices. 

We end by showing how we can deduce inversion matrices from Jacobian matri- 
ces; this was the main subject of 

Example 3.3 (Inversion matrices from jacobian matrices). In |10| . it is shown that 
one can construct a hybrid matrix whose determinant is a nonzero multiple of the 
resultant, having all rows except one of Sylvester style. In we show that the 
maximal minors of the Sylvester part of this matrix are subresultants and that we 
can solve the inversion problem by using them. 

This parameterization is extracted from [181 example Pi]; 

ti tf — ti + 1 ,2 , , 

By considering Fi := [ti + ^2)^1 ~ ^i, ^2 := + l)-'^2 - {t\ - ti + 1) and F3 := 
X3 — t\+t2, we get the following subresultant matrix: 

( X2-I 1 -1 \ 

X3 0-1-1 

Xi - 1 Xi 

Xi - 1 

\0 Xi-l Xi j 

All the maximal minors of this matrix are subresultants. Using theorem^ 
solve the inverse problem. We obtain: 

_ At, _ Xi(X2 - I-X3) _ At, _ X2X1 - X2 - Xi + 1 - XiX-i + X3 

Ai X2X1 — 1 — X2 ' Ai X2X1 — 1 — X2 

4. IMPLICITIZATION MATRICES 

We keep the notation of section |2| where we developed a matrix-based approach 
to the properness and inversion problems. In this section we introduce a new kind 
of matrices, that we will call implicitization matrices. We will prove that, when 
it exists, such a matrix characterizes the properness of the parameterization cj). 
Moreover, if is proper we can extract from it an inversion matrix of 0, as defined 
in definition 13. II 

Hereafter we will denote by F{Xi,X2, X3, X4) G K[Xi, X2, X3, X4] the implicit 
equation (which is actually defined up to multiplication by a non-zero constant in 
K) of 5 C P''. Recall that it is the homogeneous polynomial of minimal degree 



8 



LAURENT BUSE AND CARLOS D'ANDREA 



such that F{pi{t),p2{t),p3{t),p4{t)) = m ^[^1,^2,^3]; its degree is the degree of 
the surface S that we denote by deg(5'). 

Definition 4.1. A square matrix M is an imphcitization matrix of the parameter- 
ization (f) if it satisfies the three following conditions : 

(1) There exists an integer m > I and a subset V = {t^^ , . . . , t"''} of the 
monomial basis {t",a £ N'^ |a| = m} such that 



( ^ 




( P{t;X) \ 






Ah{t;X) 


V t"^ ) 




[ Mrf_i(t;X) ; 



where polynomials Mi[t; X_) are moving surfaces following (j) of bi-degree 
(m; di) and P{t;X_) is an arbitrary bi-homogeneous polynomial with positive 
degree in variables X_. Moreover there exists {t^^ , t^^ , f^^ ) g satisfying 
tit'^s = igtft and t2t'^^ = ist* (and hence tit^^ = taf^^ ). 

(2) det(M) = c.F{Xi,X2,X3,X4y where J G N \ {0} and c G K \ {0}. 

(3) If (f) is birational, i.e. S is properly parameterized by (f>, then 6=1. 

An imphcitization matrix of (j) is hence a square d x d matrix where c? > 3. 
Of course the name implicitization matrix comes from the condition (2) in this 
definition. The following theorem shows that implicitization matrices characterize 
the properness of the map <j) and moreover yield inversion matrices. 

Theoreom 4.2. Let M be an implicitization matrix of <f> and denote by M* the 
sub-matrix of M obtained by erasing its first column. Then the gcd of the maximal 
minors o/M* equals with p G N. 

Moreover the following statements are equivalent: 

(1) (j) is birational 

(2) P - 

(3) M* is an inversion matrix of S . 

Proof. First the fact that the gcd of the maximal minors of M* is a power of F 
follows immediately from the following equality in K[X] : 

(8) det(M) =F^ = CaAa 

a|— m 

where the Cq are the coefficients of the erased column of M, since F is an irreducible 
and homogeneous polynomial. 

Now suppose that (j> is birational (i.e. proper). Then we know that det(M) = cF. 
Looking at the formula © we deduce that p = since the Cq's have positive degree 
(recall that P is supposed to have positive degree in variables X). This implies that 
the rank of M* over the field K(S') is d — 1 (recall that d is the size of the square 
matrix M), that is to say that M* is an inversion matrix. 

Conversely, assume that p = 0, that is the gcd of the maximal minors of M* is a 
constant. Then the rank of M* over K(5) is d — 1 and M* is an inversion matrix. 
By theorem 13. 21 (f) is then birational. □ 

Corollary 4.3. // M is an implicitization matrix and </> is not birational, then 
det(M) = F\ with S > 1. 
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Proof. If (5 = 1 then p = in thcorem |4.2l and we would have that (p is birational. □ 

We now give a (non-exhaustive) hst of known imphcitization matrices; they can 
be divided into two distinct groups: the moving surfaces matrices and the resultant 
matrices. 

Moving surfaces matrices. All matrices coming from the moving surfaces method 
introduced by Sederberg |2()j can be used; they are by definition imphcitization 
matrices. A lot of recent works have extended the foundational work of Seder- 
berg. At this time, algorithms to construct an imphcitization matrix of a given 
parameterization cj) are available if: 

• (j> has no base points over or x P""^ (see [H]), 

• (j> has l.c.i. base points (plus some other technical conditions) over P^ (see 
El) or Pi X Pi (see jT]), 

• ip has no base points on a certain projective toric variety (see 

This method of moving surfaces is still under development (see for instance |Sll23)) 
and the list above will probably be extended in a near future. 
Resultant matrices. The imphcitization problem can be solved using some resultant 
computations, as it is illustrated in |^ The computation of a resultant often 
involves the construction of a matrix which sometimes is an imphcitization matrix. 
At their most typical, resultant matrix give universal formulas for a particular 
class of parameterizations. They are thus very interesting since they allow the 
design of a pre-computed inversion formula of a specified class of surfaces. Given a 
parameterization 0, resultant matrices which are imphcitization matrices of (j) are 
known if 

• if parameterizes a Stciner surface (see JHl; and also 0) 

• if has no base point on P^ x (see JHl and 22|^) 

• if has no base point in a certain toric variety (see 

In the list above we only mentioned general constructions and it should be pointed 
out that a lot of constructions adapted to particular cases exist. 

5. Some illustrative examples 

Below we exhibit three examples in order to illustrate our matrix-based approach 
to the inversion and properness problem. 

5.1. This example is taken from |2()j . Consider the following parameterization of 
a cubic surface with 6 base points: 

Pi = t?i2 + 2tl + tjts + Uit2t3 + Atlt:i + ititl + 2t2tl + 2tl, 

P2 = -tl - 2titl - 2tjt3 - tit^tz + titl - 2t2tl + 2tl, 

P3 = — i"^ — 2t^t2 — '^t\t^ — S^^^s — 3^it2^3 4" 2t^t^ — 2t\t^ — 2^2^35 

P4 = ^1 + t\t2 + t2 + ^I'-S "I" ^2^Z ~ ^1^3 ~ ^2^3 ~ ^3- 

If we pick the following moving planes: 

Ml :=tiXi+<2^2 + i3^3, 

M2 := ti{X2 + Xi) + t2{2X2 - X3) + t3(^2 + 2X4), 
M3 := ti(X3 - X2) H- t2{-Xi + 2X4) H- i3(^i - ^2), 



^see also |S] where it is explained that the rows of Dixon and Bczoutian matrices are moving 
planes, and hence that any hybrid combination of them give an implicitization matrix 
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we can construct the following matrix indexed by the mononiials ^1,^2,^3, whose 
determinant is the implicit equation of the surface: 

-Xi + 2X4 
X3 — 2X2 

If we erase the first row of this matrix, Theorem 14 . 21 tells us that the inverse of this 
parameterization equals (At^ : At^ : Atg), which itself equals 





X3 - 2X2 
X2 



-X2 — 2X4 
X; 



X2 + X4 



-X2 — 2X4 
X-i 



'X2 — X4 
Xi 



Xs ~ 2X2 
X2 



5.2. The following example appears in Consider the following parameteriza- 
tion: 

Pi =4+ htl - t2tl + ht2t3 - tlt2 - htl, 
P2 = ^3 + ^1^3 ^ ^2^3 ^ tli2tz + t\t2 ^ ^1^2i 



tit^ + ^2^3 



to 



P3 

Pi ~ ^3 

There are two moving planes and one moving quadric of degree one that follow the 
surface: 



^2^3 + ^li2^3 ^ ^1^2 + il^2- 



Ml 

M2 
M3 



h{Xi-X3)+t3iXi-X2), 

tiiX2 -X3 + 2X4) + t2{X2 + X3) + t3(-^2 
ti {X1X2 + X1X3) + t2 {X1X3 — X1X4 + X2 - 



- X3 + 2X4) 
X2X4) + 



+t3i-2X^ + Xl + X2X4 - X3X4 + Xl). 
In this case, the transpose of the following matrix is an implicitization matrix 
/ X4 — Xz X\ — X2 



X2 — X-i + 2X4 



X2 + X3 



—X2 — Xz + 2X4 



X\X2 + X\Xz X\Xz — X\X4 + X| + X2X4 —2X\ + X2 + X2X4 — X3X4 + X\ 1 



By erasing the last row, we obtain an inversion (At^ : Aj^ : A^g) where 
A*, 





X2 + X3 



Xl — X2 
-X2 — X3 + 2X4 



At, 



A, 



X3 — X4 
—X2 + X3 — 2X4 

X4 — X3 
X2 — X3 + 2X4 



Xl — X2 
'X2 — X3 + 2X4 



^2+^3 



5.3. We will compute the inverse of the surface parameterized by 

*? V 2 ^ _ ti+t2 



Xl 



Xo 



t'i+ti+i 



Note that here we are working with "affine" variables, i.e. we set X4 = ^3 = 1. We 
will use the Dixon formulation for the resultant. In order to do so, consider the 
polynomials 

Fi = {tl+tl + l)Xi-2ti, F2 = itl+tl + l)X2-2t2, F3 = itl+tl + l)X3-{tl+tl-l). 
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The Dixon matrix of the resultant of Fi , F2 , F3 is the 6x6 matrix 


















X, \ 




X2-2 








X2 





X2 




X3 


-1 


-1 


X3 





X3 







-X2 -2x1 + 2 


2X1 





-2X3 







-2X1 + 2-X2 





-2X3 





X2 





V 


2X1 


-2X3 








X2 


1 



whose colmTins are indexed by the monomials 1, ii, i2, ^i, ^1^2, ^2 (i'^ that order). 
So, it turns out that the transpose of D is an implicitization matrix, and we can 
recover the inverse by deleting any row and considering the ratios (Ai : Z 
For instance, by deleting the last row we get 





( ^1 







Xi-l 





Xi 


\ 






X2-2 







X2 





X2 






det 


X3 




-1 


X3 





X3 













2X1 





-2X3 











\ -X2 - 2X1 - 


f2 


-2X3 





X2 





/ 






( ° 












Xi 


\ 


•1 












X2 





X2 






dot 


-1 




-1 


X3 





X3 








-X2 - 2X1 - 


f2 


2X1 





^2X3 











\ 




-2X3 





X2 





1 






/ 


Xi 







Xi - 


1 





Xi \ 






X2-2 







X2 









X2 


det 




X3 




-1 


X3 









X3 









-X2 


- 2X1 


+ 2 




-2X3 









^ -X2-2X1+2 












X2 




) 




/ 







Xi-l 





Xi 
















X2 





X2 






dot 


-1 




-1 


X3 





X3 








-X2 - 2X1 - 


f2 


2X1 





-2X3 











\ 




-2X3 





X2 





i 





that is to say 



-2X3 (4X 



2X2) 



-X2{X, 



and t-j 



AX2X1 - 2X2 - 4X|) 
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